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We present a detailed theoretical characterization of the two fundamental collective resonances 
underlying the xenon giant dipole resonance (GDR). This is achieved consistently by two comple¬ 
mentary methods implemented within the framework of the configuration-interaction singles (CIS) 
theory. The first method accesses the resonance states by diagonalizing the many-electron Hamil¬ 
tonian using the smooth exterior complex scaling technique. The second method involves a new 
application of the Gabor analysis to wave-packet dynamics. We identify one resonance at an excita¬ 
tion energy of 74 eV with a lifetime of 27 as, and the second at 107 eV with a lifetime of 11 as. Our 
work provides a deeper understanding of the nature of the resonances associated with the GDR: a 
group of close-lying intrachannel resonances splits into two far-separated resonances through inter¬ 
channel couplings involving the 4d electrons. The CIS approach allows a transparent interpretation 
of the two resonances as new collective modes. Due to the strong entanglement between the ex¬ 
cited electron and the ionic core, the resonance wave functions are not dominated by any single 
particle-hole state. This gives rise to plasma-like collective oscillations of the 4d shell as a whole. 

PACS numbers: 31.15.ag, 32.80.Aa, 31.15.vj, 32.80.Fb 


I. INTRODUCTION 

The atomic xenon giant dipole resonance (GDR) has 
attracted much research interest since its discovery in 
1964 [1, 2], for it is one of the most prominent cases in 
atomic physics where many-body correlations play a con¬ 
spicuous role. The GDR appears in the photoabsorption 
cross section of xenon as a pronounced and nearly sym¬ 
metric hump centered around 100 eV, with a width of 
about 40 eV. The GDR lies in the electronic continuum 
above the 4 d ionization threshold. 

While the occurrence of the xenon GDR can be quali¬ 
tatively explained by the independent-particle model [3] , 
where a centrifugal barrier suppresses the 4 d —> ef tran¬ 
sitions near the 4 d threshold, satisfactory agreement with 
experiments requires inclusion of many-body correlations 
beyond the mean-field level [4-8]. Nowadays it is com¬ 
monly accepted that the xenon GDR must be described 
as the result of the collective excitations of at least all the 
4 d electrons, forming short-lived plasma-like cooperative 
oscillations [9, 10]. Because the GDR is a property of 
the inner-shell electrons, it is found in other atoms close 
to xenon in the periodic table and survives in molecules 
and solids [9, 10]. Similar giant resonances also prevail 
in nuclei, metallic clusters, fullcrences, etc [9, 10]. 

A considerable number of measurements have been 
performed for a precise characterization of the xenon 
photoabsorption spectrum in the XUV with perturba¬ 
tive light sources [11-14]. However, with the birth of 
various new source technologies, the old spectroscopic 
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feature of the xenon GDR continues to enthrall start-of- 
the-art experiments. For example, high-harmonic gen¬ 
eration (HHG) spectra of xenon driven by an intense 
mid-infrared laser display a striking enhancement in the 
plateau [15], which reflects the partial cross section of 
the 5 p valence shell strongly modified by the GDR [16]. 
Also, the GDR lies at the heart of the behavior of xenon 
exposed to free-electron lasers with ultrahigh XUV irra- 
diance [17-19]. Hence, it is important to fully understand 
the nature of the xenon GDR. 

Following the earliest independent-particle model [3], 
various advanced many-body theories have succeeded in 
reproducing the experimental cross section associated 
with the xenon GDR remarkably well [8, 20-23]. Nev¬ 
ertheless, a fundamental question frequently overlooked 
is what exactly are the basic collective modes that give 
rise to the spectral properties of the xenon GDR. A work 
by Wendin in 1971 [24] (see Ref. [8] for details) using 
the random-phase approximation with exchange (RPAE) 
identifies two collective resonances in this energy range. 
Another calculation by Lundqvist in 1980 [25] utilizing a 
hydrodynamic treatment of electron density oscillations 
also finds two collective modes, but one of them sits at 
an energy incompatible with experimental observations. 
In addition to the very limited theoretical predictions of 
the resonance positions, neither of these studies explicitly 
specify the resonance widths. Consequently, the nature 
of the inherent collective resonances hidden in the broad 
spectral blur of the xenon GDR still remains an unsolved 
question. 

The purpose of this paper is to provide a thorough 
characterization of the resonance substructures under¬ 
lying the xenon GDR within the framework of the 
configuration-interaction singles (CIS) approach [26, 27], 
an ab initio theory that can capture essential many- 
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body effects in light-matter interactions [28] including 
the xenon GDR [29]. We resolve two collective dipolar 
resonances residing in the spectral range of the GDR, 
with one position differing from that given by Wendin 
[8] by 15 eV. Whereas Wendin resorted to an approxi¬ 
mate condition only applicable to weakly damped plasma 
[30] to estimate the positions of the collective excita¬ 
tions, this work provides the first quantitative results 
for the resonance positions and lifetimes. In contrast to 
the conventional view that many-body correlations only 
quantitatively shift and flatten the resonance as seen in 
the photoabsorption spectra [3-7, 23] , we clearly demon¬ 
strate that many-body correlations qualitatively change 
the nature of the xenon GDR: a group of intrachannel res¬ 
onances splits into two far-separated resonances as soon 
as we switch on interchannel interactions involving the 4 d 
electrons. Since the resonance lifetimes are very short, 
the resonances strongly overlap and appear as one big 
hump in the photoabsroption cross section [8]. In con¬ 
trast to the plasma-type treatments used in Refs. [8, 25], 
the full many-body wave functions are directly obtained 
through our CIS approach. As the wave functions of 
the two resonances cannot be expressed by any single 
particle-hole state, we concretely show that they are in¬ 
deed new collective modes [8-10]. 

In this work, the isolation of the resonance substruc¬ 
tures is consistently accomplished by means of two com¬ 
plementary and general methods implemented using CIS. 
The first, time-independent approach provides a compre¬ 
hensive characterization of all the resonance properties 
by directly diagonalizing the many-body Hamiltonian us¬ 
ing the smooth exterior complex scaling (SES) technique 
[31-33]. Complex scaling [34] has been used to solve the 
electronic resonance problem for few-electron atoms [35- 
38] and molecules [39]. However, it has not yet been used 
to address collective resonances in many-electron atoms. 
The second, time-dependent approach involves a new ap¬ 
plication of the time-frequency Gabor analysis [40, 41] to 
the autocorrelation function of a wave packet. It is a 
common routine in molecular dynamics to look for reso¬ 
nance energies in the Fourier domain [42-44]. Nonethe¬ 
less, our analysis in the combined time-frequency domain 
not only shows improved performance in disentangling 
strongly overlapping resonances, but also supplies an ap¬ 
pealingly intuitive view on the time evolutions of various 
wave-packet components. 

The remainder of this article is structured as follows: 
Sec. II presents the theoretical tools. Sec. IIA first lays 
the foundation of our many-body CIS scheme. Sec. IIB 
and Sec. IIC explain, respectively, the time-independent 
SES and time-dependent Gabor procedures to access 
multiple resonances. In Sec. Ill we apply the methods 
to the xenon GDR, with the computational details in 
Sec. Ill A. The results of the SES and Gabor approaches 
are discussed separately in Secs. IIIB and IIIC. Restric¬ 
tions imposed on the electronic-configuration space in 
Secs. IIIB and IIIC are justified in Sec. HID. Finally, 
Sec. IV concludes the study with a future outlook. Fur¬ 


ther numerical evidence indicating the consequence of 
finding resonance poles with the approximate condition 
used by Wendin [8] is provided in the Appendix. 

Atomic units (a.u.) are used throughout the paper 
(|e| = m e = h = 47reo = 1) unless otherwise stated. 


II. THEORETICAL METHODS 
A. CIS theory 

In this work, the many-electron Schrodinger equa¬ 
tion is treated within the CIS framework, an ab initio 
theory that allows one to encapsulate essential many- 
body physics beyond the mean-field Hartree-Fock picture 
[28, 45]. Our implementation of the CIS method has been 
successfully applied to a wealth of physical phenomena 
of many-electron atomic systems interacting with light 
fields [28], including perturbative [29, 46, 47] and nonper- 
turbative [16, 48, 49] multiphoton processes with photon 
energies from the x-ray regime down to the near-infrared 
regime. Particularly, the ability of CIS to reproduce im¬ 
portant features of the experimentally observed xenon 
GDR is demonstrated in Ref. [29]. In the following, we 
outline the formulation of our CIS approach. Further 
details can be found in previous publications [26, 27, 50]. 

The nonrelativistic Hamiltonian for an IV-electron 
atom in the absence of external fields can be generally 
written as 



n^n 


=■ H 0 + Hi (1) 

where p„ and r„ are the momentum and coordinate op¬ 
erators for individual electrons, Z is the nuclear charge, 
and V MF is the mean-field potential contributing to the 
standard Fock operator [45]. The Hartree-Fock ground 
state energy E^f F is introduced to shift the entire energy 
spectrum for cosmetic purposes. The total Hamiltonian 
is divided such that Hq is merely a one-body operator and 
that all the residual two-body electron-electron Coulomb 
interactions beyond the description of the mean-field po¬ 
tential are contained in H\. 

The IV-electron Hamiltonian is represented in the IV- 
electron CIS configuration space: 

Vcis= : {K F >, ]$■>}, (2) 

which gives an ansatz for an IV-electron wave function: 

i*>=aoi$n+E <[*“>• 

i,a 


( 3 ) 
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Thus, the Hilbert space is truncated and only consists 
of the Hartree-Fock ground state |<l>(j' IF ) plus its singly 
excited configurations |$“) = c\ci | ( £ > o IF ), with c, anni¬ 
hilating an electron from an initially occupied orbital i 
and cl putting it into an initially unoccupied orbital a 
[45] . The range of the index i selects the active occupied 
orbitals from which an electron can be excited or ionized, 
i.e. the accessible channels [51], thus enabling one to test 
the multichannel character of the overall physical process 
[16, 47, 50], 

The matrix of the IV-electron Hamiltonian is then 
either diagonalized (Sec. IIB) or used in the time- 
dependent Schrodinger equation (Sec. IIC). In CIS, the 
only matrix elements that can lead to two-body effects 
are |$*). Specifically, it is the type of matrix 

elements with the indices a ^ b and i ^ j, named 
interchannel-coupling terms [51], that permits the simul¬ 
taneous change of the state of the excited electron and 
that of the ionic core, i.e., interchannel coupling leads to 
the formation of a correlated particle-hole pair [46]. Nu¬ 
merically, we can tailor the two-body nature of H\ and 
study its influences by enforcing all interchannel-coupling 
matrix elements to be zero and considering only the ma¬ 
trix elements |'I>J) with i = j. In this scenario, 

called intrachannel-coupling model [51], H\ effectively 
acts as a one-body operator: once the electron is excited, 
it can sense the potential produced by the parent ion but 
cannot modify the ionic state, which is therefore forbid¬ 
den to partake in many-body correlations [16, 29, 46, 47]. 


Solving the eigenvalue problem of the complex-scaled 
Hamiltonian with the basis set Vcis j the resonance states 
can be uniquely identified as the exposed poles situated 
above the rotated continua in the complex-valued energy 
spectrum [34, 54, 55] . It is then straightforward to obtain 
the complex resonance energy or the Siegert energy [34, 
52, 56] 

E n = a n — *r„/2, s n ,r„ G R + (5) 

as well as the wave function |$„) associated with the n-th 
resonance state. In Eq. (5), S n is the resonance position, 
and T n gives the inverse lifetime for the quasibound state 
to escape to the continuum. The detailed implementa¬ 
tion of SES within CIS using a generalized finite-element 
discrete variable representation (FEM-DVR) will be ad¬ 
dressed in a forthcoming publication [59]. 

The SES method shares with CAPs [56] the merit that 
it leaves the interior r < r$ untouched and does not per¬ 
turb the Hartree-Fock ground state (if ro and A are cho¬ 
sen suitably) [32, 34], At the same time, it eliminates 
many drawbacks of CAPs: no optimization with respect 
to a parameter is required in order to identify the reso¬ 
nance energies, and the whole transformation rests on the 
rigorous mathematical theory of complex scaling [32-34] . 

Notice that the complex-scaled Hamiltonian is no 
longer a Hermitian but a complex symmetric matrix [33] . 
As a result, the symmetric inner product (•], |-) must be 
used instead of the Hermitian inner product (-|,|-) to 
ensure orthogonality relations [27, 34, 56]. 


B. Time-independent approach to resonances: 
SES method 


A conventional procedure to access eigenstates is the 
direct diagonalization of a general Hamiltonian. Never¬ 
theless, being exponentially divergent in the asymptotic 
region renders resonance states, also known as Siegert 
[52] or Gamow [53] vectors, inadmissible elements of the 
Hilbert space of a Hermitian Hamiltonian. Standard 
techniques such as complex scaling [34, 54, 55] and the 
use of complex absorbing potentials (CAPs) [56-58] were 
thus developed to transform the wave function of a res¬ 
onance state into a single square-integrable function. In 
this paper, we adopt the SES method [31-34], a variant 
of the complex scaling technique. 

Our use of SES [33, 59] relies on an analytic continua¬ 
tion of the radial part of the electron coordinate into the 
complex plane r —> p{r) following the path of Moiseyev 
[31] in the form of Karlsson [32] adapted to the purely 
radial problem presented here: 


p(r) = r + ( e l6 


1 ) 


r + A In 


/1 + e( ro-r )/ A A 

V 1 + e r °/ x J 


(4) 


C. Time-dependent approach to resonances: 

Gabor analysis of autocorrelation functions 

Decoding the resonance substructures for a general 
quantum system can also be done through wave-packet 
propagation. The key physical quantity employed 
throughout our analysis is the time-dependent autocor¬ 
relation function defined as 

C(t)=:(*( 0)|tt(i)), (6) 

where |\P(0)) is an initial state and |4/(t)) is the freely 
evolved wave packet at a later time t. Note that the 
symmetric inner product is adopted here. This is be¬ 
cause in the time-dependent case we continue using SES, 
which effectively functions as a CAP and absorbs the 
outgoing flux reaching the end of the numerical grid 
[31, 32, 58]. For an initial state |'P(0)) orthogonal to 
|$^ IF ), the time evolution of l'F(t)) is governed by the 
CIS coefficients af(t) [cf. Eq. (3)]. Inserting Eq. (3) 
into the time-dependent Schrodinger equation with the 
complex-scaled Hamiltonian, one can derive and numer¬ 
ically integrate the equations of motion for af(t) [27]: 

j,b 


This path smoothly (depending on the parameter A) ro¬ 
tates the electron radial coordinate for r > ro about an 
angle 6 into the upper complex plane. 
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For a quantitative determination of the resonance en¬ 
ergies, we assume that, by proper preparation, the initial 
state is essentially composed of the resonance states of 
interest and all the contributions from the bound states 
and the continuum can be ignored. Expanding l'f'(O)) 
in terms of the orthonormal resonance wave functions 
|\H(0)) = J2 n an l^n), Eq. (6) then bears the following 
structure: 

( 8 ) 


The validity of Eq. (8) and the resonances that can be 
extracted evidently depend on the quality of |\I/(0)). How 
we prepare the wave packet ideal for studying the xenon 
GDR will be discussed in Sec. Ill C. 

A common strategy to infer Siegert energies from wave- 
packet propagation is to conduct a Fourier analysis and 
study the autocorrelation function in the frequency do¬ 
main [42-44]. Performing a one-sided Fourier transfor¬ 
mation on Eq. (8) [assuming C(t) = 0 for t < 0] yields 




f-\- oo 


v^F. 


= T. 


dt e iut C(t) 

+ i(uj — £i n ) 


-£■ + (u/ — ^n ) 2 


(9) 


i.e., a superposition of Lorentzian functions and disper¬ 
sive curves parametrized by the Siegert energies [60]. 

For a single resonance, the spectral distribution of the 
autocorrelation function reads 


C ( 1 ) (w) 



1 


2tt 


■p2 5 

^ + (w-Si ) 2 


( 10 ) 


which is a Lorentzian with a peak at Si and a width 
of Fi. If more than one resonance exists, |C(w)| 2 com¬ 
prises several Lorentzians and their interferences. Upon 
empirically specifying the number of resonance states, it 
is possible to retrieve the Siegert energies by numerically 
fitting \C{uj)\ 2 based on Eq. (9). 

Next, we extend the standard spectral analysis to a 
time-frequency analysis [40, 41, 61] of the autocorrelation 
function and examine its information content in the com¬ 
bined time-frequency domain. Applying a Gabor trans¬ 
formation [41, 62, 63] to Eq. (8), we can derive 



n 


Eq. (11) can be interpreted as gating the time-dependent 
signal by a Gaussian window function of width a cen¬ 
tered at t. Due to the finite size of the window function 
and the sudden turn-on of the autocorrelation function 
at t = 0, the analytical expression in Eq. (12) works as a 
reasonable approximation to Ct(w) when t^$> a. 


For a single resonance, the transient spectral distri¬ 
bution of the autocorrelation function at time t can be 
simply approximated by 


cf’M 


i e _rit e _ °’ 2(a;_E 


(13) 


i.e. a Gaussian with a peak at Si and a width determined 
by a. From the decay rate of the amplitude, Ti can be 
extracted. 

The advantage of the Gabor analysis over the Fourier 
analysis becomes apparent when multiple resonances 
come into play. In this situation, | C(w)\~ comprises sev¬ 
eral Gaussians and their interferences. Consider the ex¬ 
ample where there are overlapping broad resonances yet 
with different lifetimes. Compared to the static informa¬ 
tion conveyed by the Fourier spectrum, it is more likely to 
detect the resonances through the time evolution of the 
Gabor profile, where their competition causes dynamics 
in the frequency distribution. Quantification of the res¬ 
onance energies can be done by fitting |C t (w)| 2 with the 
help of Eq. (12) at different time steps. 

It is worthwhile to note that C(ui) is related to a 
measurable physical quantity—the photoabsorption cross 
section [29, 64]— although we usually consider its mod¬ 
ulus squared to reduce the number of irrelevant fitting 
parameters (e.g. the phase of a 2 ). Also, a real physi¬ 
cal observable—the dipole moment—can be used in the 
time-frequency analysis as well [65]. In this case, its 
frequency distribution is associated with the spectrum 
of electromagnetic radiation emitted by the system [66]. 
Lastly, it is tempting to point out the conceptual sim¬ 
ilarity between Gabor transformation and the spectro¬ 
gram measured in a pump-probe experiment [67-69], al¬ 
beit there is no real probe pulse involved in the current 
theory. 


III. RESULTS AND DISCUSSION 
A. Computational details 

The theoretical methods described in the previous 
Sec. II shall now be applied to the detailed resonance 
structure of the xenon GDR that is probed by linear 
spectroscopy using linearly polarized XUV light [1, 2,11— 
14]. The calculations are done using our XCID package 
[70]. A single set of numerical parameters is employed 
throughout our calculations to compare consistently the 
results obtained by the two approaches. 

Exploiting symmetries, nl± m is counted as one ioniza¬ 
tion channel [26, 50]. In the energy range of concern, it 
is adequate to assume that electron depopulations only 
happen from the 4 d, 5s, and 5 p orbitals [71]. However, 
for a meaningful comparison with the work by Wendin 
[8], the calculations presented in Secs. IIIB and IIIC are 
performed without activating the outer 5s and 5 p shells. 
As we will see in Sec. Ill D, these channels only cause mi¬ 
nor quantitative modifications. The HF orbital energies 











5 


are slightly adjusted to match the experimental binding 
energies. Also, the orbitals with an energy higher than 
15 a.u. are discarded to enhance the stability of the time 
propagation. 

The numerical box radially extends to a size of 150 a.u. 
and is discretized with 1800 nonuniformly distributed 
grid points with a mapping parameter of £ = 1 [27], 
from which we construct our FEM-DVR basis functions 
[59]. The SES parameters are chosen as r 0 = 20 a.u., 
9 = 40°, A = la.u. such that the Hartree-Fock mean-field 
potential remains unsealed and that the continuum is ro¬ 
tated enough to expose the resonances. The maximum 
orbital angular momentum is 3. For the time-dependent 
study, an initial state is propagated [27] at a time step 
of 0.05 a.u. for a duration of 500 a.u. to give a sufficient 
frequency resolution. In the Gabor transformation, we 
use cr = 2 a.u., which we select based on the spacing of 
the excitation energies for the two collective resonances. 
This choice represents an optimal trade-off between the 
spectral and temporal resolutions. 


B. Time-independent approach to Xe GDR: 
SES method 


The diagonalization of the complex-scaled many-body 
Hamiltonian is achieved numerically by the iterative 
Arnoldi algorithm implemented in the ARPACK li¬ 
brary [59, 72, 73]. An initial random vector is used to 
launch the iteration. Since we concentrate on the res¬ 
onance modes in the linear-response regime, the eigen¬ 
states shown below are required to have a minimum over¬ 
lap with the ground state through a dipole transition: 
| (T| D z \<f>g F )\ > 10~ 6 , where D z denotes the z compo¬ 
nent of the dipole operator [74] relative to the polariza¬ 
tion direction of the electric field in a measurement. 

Figure 1 shows the complex spectra of the energy 
eigenvalues for the case with only the intrachannel cou¬ 
plings and the case with both the inter- and intrachan- 
ncl couplings, i.e., the full-model calculation within CIS. 
Ideally, the energy spectra in both cases should follow 
the structure predicted by the Balslev-Combes theorem 
[34, 54, 55, 75]: the bound states remain on the real axis; 
the continuum is rotated clockwise by 29 degrees with re¬ 
spect to the 4d ionization threshold at 67.5 eV [76]; and 
the resonances are isolated above the continuum. How¬ 
ever, the use of an incomplete basis set results in numer¬ 
ical artifacts such as branching of the continuum away 
from the threshold [57] and a rotation angle deviating 
from 29 [77]. 

First we focus on the result of the intrachannel- 
coupling model in Fig. 1. In this case, each eigenstate 
possesses a unique hole index i and the contributions 
from different 4 d ionization channels can be easily set 
apart. Three resonances are found, one for each 4d± m 
channel. They lie fairly close to each other, forming a 
group of resonances around an energy with a real part 
ps 77 eV and an imaginary part ps —5.4 eV, which cor- 



FIG. 1. (color online). Complex energy spectra of the 
complex-scaled many-body Hamiltonian for the intrachannel 
(“intra”) and full CIS (“interH-intra”) models. Horizontal and 
vertical axes represent the real and imaginary parts of the en¬ 
ergy eigenvalues, respectively. Three close-lying resonances 
are found in the intrachannel calculation (indicated by hol¬ 
low polygons); two far-separated resonances (labeled Ri and 
R 2 ) are found in the full CIS calculation (indicated by solid 
dots). 


responds to a lifetime ps 60 as. The resonance positions 
and widths are detailed in Table I. 

In order to elucidate the origin of the small splitting 
among the intrachannel resonances, we perform another 
intrachannel calculation by approximating Hi with its 
monopole term [50]. When so doing, the resonances as¬ 
sociated with the 4do, 4d±i, and 4:d±2 channels have ex¬ 
actly the same resonance energy. Therefore, even without 
many-body effects, the electron excited from the differ¬ 
ent 4 d± m orbitals experiences different potentials owing 
to the shape of the non-spherical ionic core. This quali¬ 
tative effect, although small, clearly exemplifies the im¬ 
pact of the ionic structure beyond the description of a 
simple spherically symmetric potential (e.g. the Herman- 
Skillman [78] or Hartree-Fock [45] potentials) or even 
an angular-momentum-dependent pseudopotential [79] 
widely used to model multielectron atoms in perturba¬ 
tive [3, 80] or non-perturbative [81-83] light fields. 

Activating interchannel coupling, the group of intra¬ 
channel resonances splits into two resonance states dis¬ 
tantly located in the complex energy plane in Fig. 1. We 
have checked carefully that their positions do not vary 
with the scaling parameters, so they are not numerical ar¬ 
tifacts. In comparison with the intrachannel resonances, 
resonance i?i in Fig. 1 has almost the same excitation 
energy but a larger decay width; resonance R 2 is pushed 
away much further into the lower half of the complex 
energy plane. The splitting of the resonances highlights 
that many-body correlations are not just required for a 
quantitative agreement between theory and experiment 
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TABLE I. Siegert energies of the resonances in the intrachannel model and the full CIS model. Results are listed for the SES 
method and for the Gabor analysis. For comparison, the predictions given by Wendin [8] are included. 



SES 

a 

Gabor 

b 

Wendin [8] 

E n (eV) 

r„ (eV) 

(eV) 

r„ (eV) 

E n (eV) 

T„ (eV) 

Intrachannel: 4do 

76.3 

8.3 

76.5 ± 0.3 c 

8.2 ±0.4 c 

- 

- 

Intrachannel: 4d±i 

77.6 

13.8 

77.9 ± 0.7 c 

13.5 ±0.4 c 

- 

- 

Intrachannel: 4d± 2 

77.2 

10.6 

77.4 ±0.3 c 

10.6 ±0.1 c 

- 

- 

Full CIS: Ai 

74.3 

24.6 

80.4 ±0.7 d , 73 ±2 e 

32 ± 1 d , 17.8 ±0.4 e 

74.3 

- 

Full CIS: R 2 

107.2 

59.9 

112 ± 1 d 

47 ±9 d 

92.3 

- 


a All SES values have an error bar of 0.1 eV. This is calculated by varying over a reasonable range the numerical parameters such as the 
number of radial grid points, the maximum radial coordinate, and the SES parameters 6 and A. 
b In each calculation, the Gabor spectra are fitted numerically in a time interval if) at a time step of 3 as. This gives the fitting 
parameters S n , T n as a function of time. The error bars are then defined as the standard deviations of S n , T n over the time sequence. 
c Using Eq. (13) for one resonance in [/,,, tf] = [145 as, 363 as]. 
d Using Eq. (16) for two resonances in [/,,, tt] = [145 as, 169 as]. 
e Using Eq. (13) for one resonance in [/:,;, tf\ = [242 as, 363 as]. 


[3-7, 23], but in fact give rise to fundamentally differ¬ 
ent resonance substructures underlying the GDR. Since 
Ref. [8] does not show calculations without many-body 
correlations, our study is the first to reveal the emer¬ 
gence of the collective resonances in the GDR from the 
intrachannel resonances. 

With the interchannel interactions, each resonance 
cannot be attributed to a single ionization pathway. For 
both Ri and i? 2 , the 4 c?q 1 , and hole popu¬ 

lations [27] have a rough ratio of 1 : 2 : 1, which can 
be explained by an angular momentum analysis. Be¬ 
cause the interchannel interactions strongly couple the 
various 4d^^ hole states, it is crucial to consider the 
orbitals in addition to the one aligned along the polar¬ 
ization axis (i.e. 4 d ^ 1 ) for the physical processes involv¬ 
ing the GDR, e.g. the giant enhancement in the HHG 
spectrum of xenon [16]. We also compute the angular 
momentum composition of the excited electron, which 
shows a prominent /-wave character for both resonances. 
This is true in our intrachanncl calculation, too. Indeed, 
the xenon GDR is dominated by the 4 d —»■ ef transitions 
with roots in the independent-particle picture [3, 4], 

Our CIS approach gives the total many-body reso¬ 
nance wave functions, which are not attainable using the 
plasma-type treatments including RPAE [5, 24, 25]. We 
analyze |<f>i) and |$ 2 ) by decomposing them in terms of 
the orthonormal intrachannel basis set: 

|$n) = a ±rn\$4d ±m ) + | A$„), n= 1,2. (14) 

The first term in Eq. (14) contains the projections onto 
the intrachannel resonance functions |$4<i ±m ), and the 
second term symbolizes the remaining part with respect 
to the other intrachannel states. The complex weights 
a± m defined through the symmetric inner product are 
listed in Table II. For both Ri and R 2 , a§, aj_ 1 , and a± 2 
have the same order of magnitude. For R \ 1 the intra¬ 
channel resonances account for a weight of Y^m a m ~ 0-7 
out of a total norm of 1; for R 2 , they contribute a weight 


TABLE II. The complex weights a± m with respect to the 
intrachannel resonance states |4>4 d ±m ) for the collective reso¬ 
nances I'Ll) and |$ 2 ) in the SES calculations. 



Ri 

R 2 

a\ 

0.175 - 0.033* 

0.075 ± 0.042* 

a± 1 

0.345 ± 0.009i 

0.106 ±0.003* 

a± 2 

0.198-0.018* 

0.081 ±0.028* 

a m 

0.719 - 0.042* 

0.262 ± 0.073* 


rs 0.3. This means that the interchannel interactions do 
not only mix all the intrachannel resonance states, but 
also mix in continuum states to form the new resonances. 
At this stage, we clearly see in our Cl language why we 
can refer to R± and R2 as new collective dipolar modes: 
they are entangled particle-hole states involving the var¬ 
ious 4d±^ hole states and do not resemble any single in¬ 
trachannel resonance wave function. In RPAE, collective 
excitations are not directly defined by the many-body 
wave functions themselves, but by a coherent sum over 
different particle-hole states in evaluating the dipole [9] 
or dielectric function [8, 30] matrix elements. Note that 
both collective states have a significant overlap with the 
Hartree-Fock ground state via a dipole transition, with 
($i| D z « 1.1 and ($ 2 1 D z \$o F ) ~ T9 for R 1 and 
i? 2 , respectively. 

The positions and widths of the resonances in the full 
model can be found in Table I. The resonance positions 
calculated by Wendin [8] are also listed for comparison. 
The resonance position of R\ agrees perfectly with that 
given by Wendin. The position for R 2 differs from his 
number by 15 eV, but both positions are compatible with 
the spectral blur observed in the experimental cross sec¬ 
tion. 

The most likely reason for the discrepancy of E 2 pre¬ 
dicted by Wendin and our result is the approximate con¬ 
dition Wendin used to find collective excitations from his 
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FIG. 2. (color online). Autocorrelation functions C(t) as a 
function of time for (a) the intrachannel model (“intra”) and 
(b) the full CIS model (“inter+intra”). Solid curves for the 
real part and dashed curves for the imaginary part of C(t). 

effective dielectric function. In principle, a collective res¬ 
onance corresponds to a complex frequency where both 
the real and imaginary parts of the many-body dielectric 
function simultaneously become zero [30, 84]. An esti¬ 
mated resonance position can be found by determining 
along the real energy axis a root for the real part of the 
dielectric function, provided the damping of the true reso¬ 
nance is sufficiently small [30]. As one is dealing with two 
rather broad resonances in the case of the xenon GDR 
(particularly for R 2 ), this approximate condition, which 
is adopted by Wendin [8], is not strictly applicable. In 
the Appendix, we demonstrate how this simplified con¬ 
dition of finding the zeros of the dielectric function can 
result in Siegert energies that deviate substantially from 
the true resonance poles. Based on this argument, the 
Siegert energies provided by the present study are most 
likely to be more reliable. 


C. Time-dependent approach to Xe GDR: 

Gabor analysis of autocorrelation functions 

To probe the resonances associated with one-photon 
absorption, the initial wave packet can be conveniently 
set as 

|T(0)) = D z |< F ) . (15) 

This is equivalent to creating a wave packet via a delta- 
kick pulse polarized along the 2 axis. Since it is well 
known that the xenon GDR exhausts all the oscillator 
strength in the XUV [9, 10], the wave packet prepared in 
this way is largely composed of the relevant resonances, 
and its autocorrelation function C(t) is expected to take 
the form assumed in Eq. (8). 

The wave packet subsequently undergoes field-free re¬ 
laxation. Fig. 2 plots the time evolution of the complex- 
valued auto-correlation function for both the intrachan¬ 
nel and full CIS models. The raw data in both cases 
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co(eV) 


FIG. 3. (color online). Autocorrelation functions |C(cu)| 2 in 
frequency domain for (a) the intrachannel model (“intra”) 
and (b) the full CIS model (“inter+intra”). Solid curves for 
the data and dashed curves for the fit with Eq. (10). 


look like a simple damped oscillator without much struc¬ 
ture apart from some spikes in the beginning. This sug¬ 
gests that there are some dynamics that rapidly disap¬ 
pear. Including the interchannel interactions causes the 
autocorrelation function to attenuate faster and to ring 
at a higher frequency. 

The above features are more pronounced if we look at 
the auto-correlation functions in the frequency domain as 
illustrated in Fig. 3 [85]. The Fourier transform |C'(co’)| 2 
in each calculation shows one smooth peak, accompa¬ 
nied by Rydberg series preceding the 4 d ionization en¬ 
ergy [29]. The linewidth of the Rydberg states is narrow, 
since a 4hole decays on the femtosecond time scale 
[86]. Although the SES method yields three intrachannel 
resonances, they cannot be distinguished and really act 
as a group here. Switching on the interchannel couplings 
broadens and weakens the peak as well as displaces it 
to a higher frequency, similar to what is seen in the pho¬ 
toabsorption spectra [29] . In the full model, the two reso¬ 
nances given by the SES analysis also cannot be resolved 
in the Fourier domain. 

We can extract the effective Siegert energy for the sin¬ 
gle peak in Figs. 3(a) and 3(b). The data are fitted 
numerically with Eq. (10) utilizing the nonlinear least- 
squares Marquardt-Levenberg algorithm. This yields 
(H,r) = (79.0 eV, 12.9 eV) for the intrachannel model 
and (96.3 eV, 38.9 eV) for the full CIS model. Notice 
that |C(w)| 2 in the full CIS model is relatively poorly 
described by its Lorentzian fit and is more asymmetric, a 
hint to the multiple resonances behind the huge spectral 
hump of the xenon GDR. 

Now, we are in a position to go beyond the standard 
spectral method and to investigate the xenon GDR in 
the combined time-frequency domain. Fig. 4 depicts the 
Gabor transform \Ct(co)\ 2 for both the intrachannel and 
full CIS models [87]. In addition to the information that 
is already revealed by C(t) and IC^w)! 2 , one salient new 
feature emerges: the spectral distribution in the intra- 
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FIG. 4. (color online). Autocorrelation functions |Ct(u;)| 2 in 
the combined time-frequency domain for (a) the intrachannel 
model (“intra”) and (b) the full CIS model (“inter+intra”). 


channel case dies out almost symmetrically over time, 
whereas the spectral distribution in the full CIS model 
decays asymmetrically, with the maximum shifting to a 
lower energy. In the vicinity of 150 as, one can clearly rec¬ 
ognize two frequency components in Fig. 4(b), and can 
deduce that the higher-energy one has a shorter lifetime. 

Figure 5 presents snapshots of the frequency distribu¬ 
tions |Ct(w)| 2 in Fig. 4 at consecutive time steps, where 
the characteristics we allude to can be even better vi¬ 
sualized. In the intrachannel model, Figs. 5(a) (b)(c) ex¬ 
hibit one single, decaying peak around 80eV, which arises 
from the group of three intrachannel resonances in the 
SES calculation. Upon closer examination, we find that 
the peak position gradually moves to a slightly lower fre¬ 
quency. This is in accordance with the SES study that 
the lowest-lying intrachannal resonance has the smallest 
decay width (see Table I). 

In the full CIS model, the time evolution of the tran¬ 
sient spectral distribution is fundamentally different. In 
Fig. 5(d), the initial spectrum displays a broad and 
nearly symmetric peak located around 90 eV. However, 
the spectrum soon becomes highly asymmetric with the 
maximum shifting to a lower energy, and decays faster 
than the intrachannel case. The successive dynamics 
in Fig. 5(e) vividly picture how the substructures in 
the GDR develop—two peaks can be identified, and the 
higher-lying one fades away much quicker. This is fol¬ 
lowed by Fig. 5(f), where the higher-lying mode has com¬ 
pletely vanished and only the lower-lying one remains, 
with a position similar to that of the intrachanncl res¬ 
onance. Just based on simple observations, the Gabor 
analysis intuitively illuminates how the interchannel in- 
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FIG. 5. (color online). Decay dynamics of the Gabor spec¬ 
tra \Ct(oo)\ 2 at consecutive time intervals. Panels (a)(b)(c) 
in the the left column for the intrachannel model (“intra”); 
panels (d)(e)(f) in the right column for the full CIS model 
(“inter+intra”). In each panel, the time step between two 
neighboring lines is about 3 as. 


teractions result in the damping and fragmentation of 
the resonances, as well as a rough idea of the resonance 
positions and widths. Benefitting from the fact that dif¬ 
ferent spectral components have different lifetimes, the 
Gabor analysis successfully disentangles the two funda¬ 
mental collective modes that cannot be separated by the 
Fourier analysis. 

Next, the Siegert energies are quantitatively deter¬ 
mined following the logic presented in Sec. IIC. Based 
on the a priori input from the SES study, we perform 
the analysis in the intrachannel case for each ^d± m ion¬ 
ization channel. One peak in the Gabor profiles at 
subsequent time steps is then fitted numerically with 
Eq. (13). The outcomes are tabulated in Table I, and 
are in excellent agreement with the SES results within 
the error bars. Particularly, the Gabor analysis cap¬ 
tures the minute splitting trend of the resonance energies, 
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i.e. ^4d ± i > ^4d ±2 > “4c?o an d r 4 d ± i > r 4d±2 > r 4do . 

For the full CIS model, the fitting process is dissected 
into two stages. At the first stage, roughly corresponding 
to the time interval shown in Fig. 5(e), two resonances are 
singled out. Resorting to Eq. (12) with n = 2, the Gabor 
spectrum has the approximate analytical expression: 


cf>M 


~ /mM +/2,t(w) 

+ 2 [/i,t(w)/ 2 ,t(w)] 3 /i_2,t(w), (16) 


where 


, 2 | 2^ ? _ rnte _ ff2(w _ En)2 


fn,t(u) := a n e 4 e " e " 


/i- 2 ,t(w) :=cos 


tr 2 (r 2 - ri)o 


+ (r 2 - rot 


i (r 2 s 2 -r 1 s 1 ) 


+ (<(>2 - 0l) 


■ (17) 


In Eq. (16), the first two terms are the individual contri¬ 
butions from R\ and R 2 , and the third one is their in¬ 
terference. In Eq. (17), 0 n denotes the phase of a 2 . The 
data are fitted with the above formulae, and the Siegert 
energies are given in Table I. The error bars here are big¬ 
ger than those in the intrachannel case, and the extracted 
Siegert energies deviate from the SES ones, especially for 
the decay width of the faster-decaying R 2 . This is possi¬ 
bly due to the increasing difficulties in the nonlinear fit¬ 
ting procedure (particularly from the interference term). 
Also, in order to arrive at the analytical expression for 
|Ct(w)| 2 , Eq. (8) assumes that all the contributions to 
C{t) from the bound states and the continuum can be 
neglected, which works worse if the resonances decay rel¬ 
atively fast. 

At the second stage, which nearly coincides with the 
time interval shown in Fig. 5(f), only one resonance is 
seen. Using Eq. (13), we produce another Siegert en¬ 
ergy for i?i in Table I. It does not fully agree with that 
retrieved at the former stage, especially in terms of the 
resonance width. The two-resonance model used at the 
first stage is only applicable in a short period of time, 
where the resonance parameters certainly cannot fluctu¬ 
ate too much. Hence, this discrepancy reflects further 
numerical instability in the fitting parameters that can¬ 
not be entirely represented by the previously calculated 
uncertainties. The Siegert energy acquired at this sec¬ 
ond stage seems closer to the SES one. Nevertheless, 
as the Gabor spectrum in this time interval has a fairly 
weak amplitude, the contribution from the Rydberg se¬ 
ries (after the filtering) inevitably kicks in, which lowers 
the effective energy position and width for R\. 

The most appealing feature of the Gabor analysis 
is that it provides an intuitive dynamical view on the 
competition between various spectral components, which 
may be connected to what is measured in a pump-probe 
experiment [67-69]. However, it is apparent that the 
Gabor method cannot quantify Siegert energies as accu¬ 
rately as the SES approach. Considering the deviations 


from the SES results, the energy uncertainties, and the 
discrepancy between the resonance energies for R\ ex¬ 
tracted at two different stages, the Gabor analysis gives 
an overall energy resolution ss 10 eV. 


D. Influence of 5s and 5p orbitals 

In the above Secs. Ill B and III C, the active ionization 
channels lie in the 4 d shell; the outer 5s and 5 p shells are 
frozen. This is an assumption made in Ref. [8] as well. 

Our SES calculations show that including the 5s and 
5 p± m channels leads to no qualitative but only minor 
quantitative modifications to the previous discussion. 
Hence, the xenon GDR mainly stems from the many- 
body correlations involving the ten 4 d electrons [9], and 
the 5s and 5 p electrons are only small admixtures. The 
Siegert energies for the three intrachannel resonances re¬ 
main the same as in Table I. In the full CIS model 
with active 4d, 5s and 5 p shells, the Siegert energies 
are slightly revised to (Si,Ti) = (73.4 eV, 24.7 eV) and 
(S 2 , r 2 ) = (111.8eV, 58.2eV). Note that the higher-lying 
R 2 is more sensitive to the effects of the outer shells. 


IV. CONCLUSION 


In this paper, we disentangle two fundamental collec¬ 
tive dipolar resonances that cannot be resolved in the 
photoabsorption cross sections associated with the xenon 
GDR. In extension of Wendin’s pioneering work [24] , we 
achieve a complete theoretical characterization of the 
resonance substructures by two complementary meth¬ 
ods within the framework of the CIS theory. It is very 
likely that the Siegert energies given by the current study 
are more accurate than those given by Wendin, as our 
methodology for finding the resonance poles is not limited 
to weakly damped oscillations. The time-independent 
SES approach demonstrated here is the first example 
of treating collective resonances in multielectron atoms 
with the complex scaling technique. The time-dependent 
Gabor analysis extends the standard Fourier analysis to 
the combined time-frequency domain, such that strongly 
overlapping resonances living on different time scales can 
be more easily separated. 

Our work provides a deeper insight into the nature 
of the GDR: the group of three close-lying intrachannel 
resonances splits into two far-separated resonances upon 
the inclusion of interchannel couplings primarily involv¬ 
ing the 4 d, electrons. The two resonances are new col¬ 
lective modes in the sense that they must be written as 
a superposition of various particle-hole wave functions. 
When the excited electron is still near the ion, a strongly 
entangled particle-hole pair is formed. This leads to the 
strong mixing of the various 4dlj] m ionic states, the entire 
4 d shell thus exhibiting collective plasma-like oscillations 
as a whole. 
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We specify the Siegert energies for the two collec¬ 
tive resonances. However, the exact values need further 
theoretical refinement. The CIS theory only contains 
one-particle-one-hole configurations (in addition to the 
Hartree-Fock ground state). Hence, real and virtual dou¬ 
ble excitations [8, 9, 51] are among the physical processes 
outside the scope of the current study. Nonetheless, since 
TDCIS (in the velocity gauge) produces a peak position 
in good agreement with the experimental photoabsorp¬ 
tion cross section [29], we expect that inclusion of dou¬ 
ble excitations would not affect the resonance parameters 
substantially. 

Finally, we note that a recent experiment at the free- 
electron laser FLASH, using an XUV nonlinear spec¬ 
troscopy technique, has provided the first direct evidence 
for the two collective dipolar resonance states associated 
with the xenon GDR [90]. Thus, it may be expected that 
experiments of this type will provide an opportunity to 
test the predictions presented in this paper. 
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Appendix: Consequence of using the approximate 
condition to find zeros of dielectric functions 


structure [84]: 


e(E) 


1 

1 — ( °i _|_ h _’ 

1 \E-E 1 ^ E-E 2 ) 


(A.l) 


where ai and a 2 are two complex numbers. 

To extract Wendin’s dielectric function for the xenon 
GDR, we fit the real part of his data with Re(e(S)) in 
Eq. (A.l). The a n are treated as fitting parameters, and 
the E n as known constants with the values given by the 
SES method in Table I (as mentioned in Sec. Ill B, we 
believe that our Siegert energies are closer to the true 



FIG. 6. (color online). Real part of the many-body di¬ 
electric function Re(e) along the real energy axis S. Solid 
curve for Wendin’s data [8] and dashed curve for the fit with 
Eq. (A.l). The two resonance positions determined in the 
spirit of Wendin’s work are labeled S( and S 2 . 


As briefly explained at the end of Sec. IIIB, the ap¬ 
proximate condition invoked by Wendin [8] to find ze¬ 
ros of the many-body dielectric function can result in 
Siegert energies that substantially deviate from the true 
resonance poles. In this section, we provide numerical 
evidence in support of the above statement. 

A collective resonance obtained from diagonalizing the 
complex-scaled many-body Hamiltonian corresponds to 
an energy E n = S n — iT n /2 in the complex energy 
plane where both the real and imaginary parts of the 
many-body dielectric function e(E) simultaneously van¬ 
ish [30, 84], In the limit of T n —> 0, this exact condition 
is reduced to finding the roots of the real part of e(3) 
along the real energy axis 3, T = 0 [30]. 

For a dielectric function with two poles in the energy 
range of interest, it must follow the simple analytical 


ones). The real part of the reconstructed dielectric func¬ 
tion (with a\ = —10.0 — 19.7* and a 2 = —1.6 + 6.4*) is 
shown in Fig. 6, and nicely describes Wendin’s result. 

The real part of e(S) we retrieve passes through the 
real energy axis at 3] = 74.5 eV and E' 2 = 91.0 eV, which 
do not coincide with the excitation energies of the true 
collective poles Si = 74.3 eV and S 2 = 107.2 eV. This 
strongly indicates that the approximate condition of find¬ 
ing the roots of Re(e(S)) as Wendin did does not suffice 
to provide accurate predictions for the short-lived collec¬ 
tive resonances. In particular, the estimated S 2 for the 
shorter-lived i? 2 amounts to an error of 16.2 eV. In our 
opinion, the most consequential approximation Wendin 
made does not lie in the way he constructed the dielec¬ 
tric function, but in the way he searched for the collective 
poles, which is suitable only for weakly damped oscilla¬ 
tions [30]. 
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